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Abstract 

We develop velocity-dependent models describing the evolution of string networks that involve 
several types of interacting strings, each with a different tension. These incorporate the formation 
of Y-type junctions with links stretching between colliding strings, while always ensuring energy 
conservation. These models can be used to describe network evolution for non-abelian strings as 
well as cosmic superstrings. The application to Zjv strings in which interactions are topologically 
constrained, demonstrates that a scaling regime is generally reached which involves a hierarchy of 
string densities with the lightest most abundant. We also study hybrid networks of cosmic super- 
strings, where energetic considerations are more important in determining interaction outcomes. 
We again find that networks tend towards scaling, with the three lightest network components being 
dominant and having comparable number densities, while the heavier string states are suppressed. 
A more quantitative analysis depends on the precise calculation of the string interaction matrix 
using the underlying string or field theory. Nevertheless, these results provide further evidence 
that the presence of junctions in a string network does not obstruct scaling. 
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I. INTRODUCTION 



Much of the interest in cosmic strings was lost with the reahsation that they cannot be 
the dominant seeds for structure formation in the Universe . However, their appearance 
in many cosmological situations forces one to consider them as subdominant contributors. 
Recent theoretical work in brane inflation and SUSY GUTs [6,], as well as the 

potential for observing strings through gravitational lensing 0, y, 3 and the CMB , have 



contri 
Refs. 



3uted 
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;o a significant revival of interest in the subject (for reviews see for example 



12 



13 



14|). More recent work [15|, [l6| suggests that CMB data allow a 10-15% 



contribution from cosmic strings, arguably favouring a scale- invariant spectrum plus strings 
over a tilted spectrum without strings. 

Particularly interesting is that cosmic strings can be produced in the final stages of 
brane inflation 



It 



18 



cosmological distances [21 



22 



20], as D-branes and fundamental strings stretching over 



231]. These objects, often referred to as 'cosmic superstrings', 
can have different phenomenology than ordinary field theory strings. This opens up the 
possibility that cosmic string observations could yield information about physics at the string 



scale, and provides new ways of constraining various brane inflation models 
Cosmic superstrings have small tensions (in the range 10^^^ < Gfi < 10~^ 
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J26|,12 



reconnect with probabilities that can be significantly less than unity 



28 



28|) and 



2^ . This leads to a 



reduced intercommutation rate, resulting in an enhancement of the predicted string number 



density tod ay 22l . l28l . 1301 ]. Cosmic superstring networks can consist of more than one type 
of string 0, 23], which can zip together to produce trilinear vertices with links (or better. 



in the terminology of Ref. 



31 



similar to non-abelian strings 



, zippers) stretching between them. In this respect they are 



31 



32 



33 



34| . In contrast to abelian string networks, which 



have been shown to always reach a scaling solution 35, 



could in principle be frustrated 



36 



37 



38l . l39l ] , non-abelian strings 



33]. 



The evolution of non-abelian string networks was further studied in 



31 



34| where some 

evidence for scaling behaviour was found. There have also been more recent attempts to 



40 



41 



42 



45|, 



46|, 



model non-abelian strings both analytically and numerically 

favouring scaling (see also Ref. {4^ where scaling was found for domain wall networks with 
junctions). In this paper we present a class of velocity-dependent models for non-abelian 



string evolution, developed under significantly different assumptions than those of Ref. 



4l|. 
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In particular, rather than associating a different energy density to each type of string, 
while keeping a single correlation length and average velocity for all types, we associate a 
different correlation length and velocity to each string type and assume a Brownian network 
structure, in which correlation lengths also quantify string energy densities. We consider 
two basic categories of interactions between strings of different types, namely the coalescence 
along their own length (zipping) that produces a 'zipper', and the creation of a new segment 
(a 'bridge') of different string type as the strings pass through one another. The relative 
importance of these inter-related dynamical mechanisms will depend on the details of the 
field theory model under study, and more general interactions can be described as relative 
combinations of these two limiting cases, as we shall comment later. Our network evolution 
models are constructed in a manner close to the traditional (Kibble) approach, where the new 
terms corresponding to the production of bridges/zippers, are introduced through energetic 
considerations of string interactions within a certain volume element. We express our models 
in terms of the 3D energy density of strings, thus ensuring the use of well established 
techniques in calculating the energy loss of the long string network. We can then directly 
compare our results to both the usual (3D) VOS model 3^ and numerical simulations 



38, 



48|. 



Central to our models is ensuring that energy is conserved by balancing the energies 
corresponding to the produced and lost string lengths and/or the kinetic energies of the 
interacting and produced string segments. In this way we avoid the possibility of finding 
'spurious scaling' by artificially throwing away energy from regions of intercommutation. 
Realistic collisions can of course be inelastic, in that some of the energy can be released in 
the form of radiation. This provides an additional energy loss mechanism, which facilitates 
the reaching of a scaling regime. In the absence of a quantitative understanding of the 
importance of such an energy loss mechanism, we take a conservative approach and assume 
energy balance at the interaction regions. We find that even in this case, where this poten- 
tially important energy loss mechanism is effectively switched off, multi-tension networks 
with junctions can still reach scaling. 

The structure of the paper is as follows: In section [ITl we briefly review some basic facts 
about cosmic string evolution both in the abelian and non-abelian case. In section IIIII we 
move on to construct a class of velocity-dependent models for non-abelian string evolution. 
Application to Zn strings is presented in section \W\ whereas in section |V] we consider the 
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case of cosmic superstrings. Our conclusions are summarised in section |VT1 



II. COSMIC STRING EVOLUTION 



A. Abelian Strings 



Most well-studied examples of cosmic strings are those arising in theories with a broken 
f/(l) symmetry, like the Nielsen-Olesen vortex lines of the abelian-Higgs model. Under- 
standing the interaction of strings when they intersect is crucial for modelling their cosmic 
evolution. For such U{1) strings, crossing has two topologically allowed outcomes: either 
the strings pass through one another, or they intercommute (reconnect) by exchanging part- 
ners. Which situation is actually realised is a dynamical question. This problem has been 



extensively studied numerically in field theory 49|, |50(], where it was found that the strings 
reconnect with probability of order unity 8^1. Some recent analytic results have also been 



presented in Ref. Sjj. Having quantitative understanding of string-string interactions is 
necessary for modelling the cosmological evolution of strings. Below we review some basic 
results. 



1. Basics 

Numerical simulations of string formation in field theory suggest that strings formed after 
cosmological phase transitions have, to a good approximation, the shape of random walks. 
This allows one to describe a string network by a characteristic length L, which determines 
both the typical radius of curvature of strings and the average interstring distance. There 
is typically one string segment of length L in each volume so that the energy density of 
the network can be defined as 

pL 

p = U = T-- (1) 

As the network evolves, strings collide with each other or curl back on themselves creating 
small loops, which oscillate and radiatively decay. Via these interactions enough energy is 
lost from the network to ensure that the strings do not dominate the energy density of the 



universe |81l] . An approximate energy loss rate equation can be written [35 1 



where the first term is due to Hubble expansion (a = a{t) being the scalefactor of the 
universe) and the second term models the energy lost through string interactions and the 
formation of loops. Such interacting networks are known to evolve towards a so-called 
'scaling' regime in which the characteristic length L stays constant relative to the the horizon 



d 



H 



t 



35| . This can be seen by setting L = 'y{t)t and substituting ([T]) into ([2]) to obtain 



7 



1 = - 2{f3-l) + - 



7 



2t 



7. 



(3) 



The expansion exponent (3 is related to the scalefactor a{t) by a{t) oc and is equal to 
1/2 and 2/3 in the radiation and matter eras respectively. Equation has indeed a scaling 
solution 

7= [2(1-/3)]-^ (4) 

demonstrating that the characteristic length asymptotically reaches a constant value with 
respect to the horizon L ~ t. If one starts with a high density of strings, intercommuting 
will produce loops reducing the energy of the network, whereas if the initial density is low 
then there will not be enough intercommuting and 7 will decrease. Given enough time, the 
two competing effects of stretching and fragmentation will always reach a steady-state and 
the scaling regime will be approached. 

The above discussion captures the key physical processes involved in string evolution, but 
for a more quantitative study one needs models of higher sophistication. String evolution has 



37|, 



38 



481], verifying 



been extensively studied numerically through Nambu-Goto simulations 
the above picture, but a number of analytic models have also been developed. These include 
a velocity-dependent one- scale (VOS) model Q 13, 



functional approach 



53l | , a 'kink-counting' model [SJ, l55| , a 



561 ]. a 'three-scale' model [57| and a 'wiggly' model 



In the following 



we will consider the VOS model, as it is the simplest of these and has been shown to be in 

□ 

very good agreement with numerical simulations [591 ]. 



2. The VOS model 



The velocity-dependent one-scale (VOS) model is a simple analytic model, depending 
on one free parameter only, which agrees quantitatively with Nambu-Goto string evolution 
simulations. From the theoretical point of view it is well-motivated and can be obtained 
directly from the Nambu-Goto action, by performing a statistical averaging procedure on 

5 



the string equations of motion and energy momentum tensor along the network [52]. In 
particular, from the energy momentum tensor one can obtain an evolution equation for 
the string energy density p, whereas the Nambu equations of motion yield a 'macroscopic' 
equation for the evolution of the typical rms velocity v of string segments. 
The resulting equations are: 

p=_2«(l + „>_£^ (5) 

Qj Li 

(6) 

The second term in equation ([5]) is a phenomenological term which takes into account 
the energy losses through creation of loops. This depends on the loop production parameter 
c, related to the integral of an appropriate loop production function over all relevant loop 
sizes Q]. In equation ([H]), -R is the average radius of curvature of strings in the network 



and k the so-called momentum parameter (see later discussion) introduced in Ref. [52] • For 
a Brownian network, and within the VOS assumptions, the average radius of curvature R 
can be taken equal to the correlation length L = jt. 

Using ([1]) and setting L = ■y(t) t as before, we obtain the following equation instead of 
©: 

2^1(2,(1 + „V 2 + ^). (7) 

This is of the same form as ([3]) but has an extra correction term (3v'^ accounting for 
redshifting of velocities due to cosmological expansion. It also includes the parameter c, 
;he value of which can be extracted from numerical simulations and it is of order unity 



32 



531]. Possibly reduced intercommuting probabilities relevant to cosmic superstrings. 



or non- intercommuting for high velocities in type II strings 60|], correspond to a smaller 
effective value of c. 

The system (jMZl) has the scaling solution 

^ 4/?(l - /5) P{k + c) ^ ' 

in terms of the expansion exponent /3, the loop production parameter c, and the momentum 
parameter k which is a measure of the smoothness of the strings. An accurate ansatz for 



the momentum parameter has been proposed in 



39] 



, , 2^2 (l-Sv^\ 



which incorporates the 'Virial' condition v"^ < 1/2, observed in Nambu-Goto simulations 
32|. 

With this ansatz the VOS model depends on one parameter only, c, which, as mentioned 
above, can readily be extracted by comparison to numerical simulations. Remarkably, by 
adjusting this parameter once only, one can closely model high resolution Nambu-Goto 
simulations throughout cosmic history. 



B. Non-Abelian Strings 

Consider a situation where a gauge symmetry group G is spontaneously broken to a 
subgroup H. This process can give rise to topological defects classified by the homotopy 
groups of the vacuum manifold G/H. In particular, string defects can occur if the vacuum 
manifold is not simply connected that is 7Ci{G/H) ^ I, where tti is the fundamental group 
and I the identity. In the example of U{1) strings the symmetry breaking is U{1) I 
so TiiiG/H) = Z, the group of integers, corresponding to the integer winding number of 
U{1) strings. For a connected and simply connected gauge group, the fundamental theorem 
implies tti^G/H) ~ -kq^H), where tcq^H) is the set of disconnected components of H. Thus, 
for a simply connected G, strings can occur if H is disconnected. 



The fiux / of a string is given by a path-ordered exponential of the gauge field Gjj 



f = Pexpl^i j^^A^dx^^ , (10) 

where 7 is a closed, oriented path, which encloses the string without passing from the string 
core and without enclosing any other strings. The requirement that the Higgs field which 
drives the symmetry breaking is invariant under parallel transport along such a closed path 
forces the fiux to be an element of the unbroken group H. If H is non-abelian then the 
above definition is not gauge invariant: under a gauge transformation hy h E H the fiux / 
changes to hfh~^. However, two paths which start and end at the same point p and which 
can be continuously deformed to each other necessarily have the same fiux. 

The non-abelian structure of H gives rise to a certain degree of ambiguity in defining the 
fiux of a particular string. Consider a string A which is entangled to another string B as in 
Fig. [TJ One can calculate the fiux a on one side of A by considering the path a, and the 
fiux on the other side a' by considering path a'. Since strings A and B are entangled, path 
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p 



FIG. 1: The paths a and ol on the two sides of the strings are not homotopically equivalent and 
as a result the corresponding fluxes are related by conjugation with /?, that is ol = f5aj3~^. 

a cannot be deformed to a' without encountering string B. Therefore, they correspond to 
different elements of the fundamental group and the fluxes a and a' are generally different. 
In fact they are related by conjugation with the flux (3 (of string B, corresponding to path 
h) that is 

a' = l3al3-\ (11) 

which gives a' = a in the abelian case. This is not merely a mathematical curiosity, but 
a physical effect that can be interpreted in terms of a long-range interaction between non- 
commuting strings, the so called 'holonomy interaction' GlJ, l62, l63|]. Thus, the flux of a 
particular string is not associated to a unique element but can be described by several 
distinct conjugate elements oi H . It follows that strings with fluxes in the same conjugacy 
class must have the same tension. 

Non-commutativity of H changes the nature of string interactions, thus giving rise to 
an interesting property of non-abelian strings: the ability to form entangled networks. In- 
deed, two strings carrying non-commuting fluxes can neither pass through each other nor 



reconnect, as both possibilities would violate flux conservation [6J, |65j. The collision of two 
such strings leads to the formation of a new string segment stretching between the two, 
which carries flux given by the commutator of the initial fluxes. A related process is that 
of branching, that is the splitting of a string with flux c = a6 to two strings with fluxes a 
and h. This is not however special to having a non-abelian unbroken group if; it can also 
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occur for strings arising through the process U{1) — > Z^, with > 3 [3J, [Gy, l67| . Also, for 
type 1 U{1) strings, coUiding segments can zip for small enough crossing angle and velocity 
of approach 



The evolution of entangled and/or branched networks is not as well understood as that 
of ordinary U{1) string networks. A certain class of non-abelian strings have been studied 



numerically in Ref. 



where the possibility of a 'frustrated' network (one that does not 



exhibit the scaling property and eventually dominates the energy density of the universe) was 



identified. Simulations of 5*3 strings were presented in Refs. [31|, l3J|, where strong evidence 
for scaling behaviour was found. Interest in the subject has recently revived in the context 
of 'cosmic superstring' networks, which also have the branching property. Recent attempts 



to model such networks both analytically and numerically can be found in Refs. 



42 



45 



40 



41 



46|. In particular Ref. [40|] suggests a way to generalise the VOS model to take into 



account the formation of links (or 'bridges') in an entangled/branched network. Ref. 41| 
proposes a further generalisation in order to model networks composed by different types of 



42, 



45 



46| take a field theory approach 



strings, each with different tension. Finally, Refs. 
to cosmic superstring models. 

The aim of this paper is to develop an analytic model for the evolution of multi-string, 
entangled/branched networks in the spirit of Refs. 40|, |4l|. Our work closely follows the 



traditional Kibble approach to modelling string-string interactions in terms of energetic 
considerations in a correlation length volume. We start developing our models in the next 
section. 



III. NON-ABELIAN VOS MODELS 

We wish to study a network of cosmic strings consisting of several different types of string, 
possibly with different tensions and non-commuting fluxes. The two mechanisms governing 
the evolution of the network are the cosmological expansion (with its associated velocity 
redshift effects) and the string-string interactions. For U(l) strings both effects are well- 
described by the VOS model, but to allow for the new possibilities arising in a rnulti-string 
network (section III Bl) this model needs to be accordingly modified. As in Ref. j4l| we aim 
to have a set of VOS equations for each type of string, with extra terms added in order to 
describe the interactions between different string types which transfer energy from one type 
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FIG. 2: The bridge-type interaction: the colUding strings form a new segment, the 'bridge', as they 
pass through each other. 



to another. 

String interactions are modelled as follows: when two strings of the same type collide 
there are two topologically allowed outcomes; either they intercommute or pass through one 
another. This will be described by the usual VOS parameter c (see Eq. (^), parametrising 
the energy lost through self-intersection of strings and the subsequent production of loops 
of the same type. However, if the collision involves different types of strings (with non- 



commuting fluxes) then both the above possibilities violate flux conservation. The co 



lision 



instead leads to the production of a new string segment, joining the two strings [82]. As 
already mentioned we choose to describe the formation of the new segment via two inter- 
related limiting mechanisms: 

(i) The strings pass through one another and in doing so they become linked by a new string 
segment of a third type (Fig. [2]). This is the so-called 'bridge' configuration in the terminol- 
ogy of Ref. [31^. In this asymptotic limit the energy to make the new link is acquired from 
the kinetic energy of the original pair. 

(a) The strings coalesce along part of their own length, in the so-called 'zipper' configura- 
tion of Fig. [31 In this limit there must be significant energy gain from the realignment and 
coalescence of the two original strings, which lose physical length to create the new segment. 

A general model will involve some combination of these two limiting mechanisms. Their 
relative importance is both a dynamical and kinematic question depending on the string 
tensions, velocities and angles between the colliding strings, as well as topological considera- 



tions. Mo del- dependent averages, such as those considered in Refs. |43|, [4^ , may be used to 
determine their relative quantitative importance. We will first consider the impact of each of 
these possibilities separately in sections llll A 21 and IIII A 3l Then we will present the general 
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FIG. 3: The zipper-type interaction: strings zip along part of their own length to form a 'zipper'. 



case with both mechanisms incorporated to describe more general string interactions. 



A. Toy Models 

Consider a simple example where we begin with a network of Brownian strings of types 1 
and 2, with tensions fii, fi2 and correlation lengths Li, L2 respectively. We can think of this 
system as two separate networks, each characterised by its own correlation length, which 
are allowed to interact with each other. Interactions between strings of the same type are 
as in the U{1) case, but when strings of type 1 and 2 collide they produce a new segment 
of type 3 with tension fi^. We assume that the energy density of segments of type 3 can be 
described by a correlation length L3, as for a Brownian network. We can then define the 
energy density ps = Hzj L\ and write an evolution equation 

P3 = -2-(l + t;3)p3 +Pl,2^3- (12) 

a 

The first two terms are as in Eq. ([5]) and describe the effects of cosmological expansion and 
loop production through self- intersection of strings of type 3. The third term models the 
production of new segments of type 3 string from interactions between strings of types 1 
and 2. 

To write down a formula for pi,2^3 we consider a correlation-length-sized segment of type 
1 string interacting with the type 2 network and take, without loss of generality, Li < L2. 
On average there is one type 2 segment of length L2 in each correlation volume Lg, so that, if 
V is the velocity of the type 1 segment relative to the (type 2) network, then the probability 
that the segment meets a string of type 2 in time 6t is: 

11 



Such a colhsion results in the production of a new type 3 segment. In analogy to the case of 
loop production, we will assume that the length distribution of the produced links is peaked 
at a value i{t). Integrating over this distribution function introduces an efficiency parameter 
d, the analogue of parameter c for the case of loop production. Considering that on average 
there are -j^j-^L^ = L\/L\ strings of type 1 in each cube of volume L\ we can write 

_ jv Li H3i{t) Ll _ dvn3i{t) 

Pl,2^3-a-j-J 73—73- t2t2 ■ 

1j2 J^2 J^2 ^1 -^1-^2 

Note that this expression is symmetric in Li, L2 as expected. 

With equation ([Hj) modelling the production of links due to interactions of strings of 
different type, and with the assumption that the produced links form a Brownian network, 
we can go on to construct simple VOS models of non-abelian string evolution. We will first 
consider toy models of networks consisting of three types of string only, which interact in a 
simple manner, and in the next sections we will progressively build more complex models. 
Before presenting the models we briefiy discuss the important issue of energy conservation. 



1. Energy Conservation Issues 

When considering macroscopic string evolution equations, it is important to make sure 
that energy is conserved by the interactions. In the 3D VOS model, the loop production 
term damps energy from the long string network, but this energy is simply transferred to the 
network of cosmic loops. However, for processes where two colliding long string segments 
produce a (long) segment of third type, it is important that the energy differences balance 
exactly, so that no net energy is produced or lost. If this condition is not met, one is in 
danger of unphysically throwing energy away, which may lead to a 'spurious' scaling regime. 
Of course, string collisions may well be inelastic, and it seems plausible that some energy 
can, indeed, be lost through particle production 0, Q or gravitational radiation. However, 
it is at present not clear what proportion of the relevant energies can indeed be damped away 
from the network through, say, massive radiation. In the case of Abelian-Higgs networks, 
there is evidence that massive radiation is subdominant compared to loop production and 
gravitational radiation decay mechanisms |69| . 

In the absence of a quantitative understanding of particle production in the present 
context, it seems dangerous to assume that the energy difference associated to the zipping, 

12 



say, process is all damped away through this decay channel. A more conservative approach 
would be to impose energy conservation during string intercommutations (in other words 
to switch off particle production), and see whether the resulting network;s can still reach 
a scaling regime. That is the approach we will follow in this paper. As we will see, for 
interactions of the bridge type, this energy balance can be accomplished by taking into 
account the slowing down of the original strings due to the tension of the bridge between 
them, or, in the zipper case, by giving kinetic energy to the produced zipper through terms 
in the corresponding velocity evolution equations. We now examine each of these situations 
separately. 



2. Model with Bridges 

Now consider the case where strings of type 1 and 2 produce type 3 segments through 
interactions of the 'bridge' type (Fig. [2]). As the strings pass through one another, a new 
segment of type 3 is produced, while the initial length of each of the original strings is 
preserved (Note that by 'length' we mean physical, not invariant, length). The tension of 
the formed link slows down the original strings, and in doing so its length increases. If i{t) 
is the average value of the lengths of the links at the end of all interactions which occurred 
at times around t, then the evolution of the energy density of the type 3 network is given 
by Eqs. f|T2|) and f|T^ . However, since the length of type 1 and 2 strings is not affected by 
the production of bridges, the corresponding evolution equations for networks 1 and 2 will 
be of the form of Eq. We will therefore have for the energy densities: 

p. = -2^{l + v^)p,-^^, ^ = 1,2 (15) 

P3 = -2-{l + V^)p3- ——+ . (16) 

3 12 

Equations (fT5|) - (|T6|) seem to violate energy conservation since for d = (flat space) and 
Cj = (no loop production), the first two equations become pi = but has a positive 
term, appearing to create energy from nothing. This energy is actually coming from the 
kinetic energies of the colliding strings (which have been slowed down due to the tension 
of the bridge stretching between them) and needs to be taken into account in the velocity 
equations. To do this, we express the string density p in terms of the proper correlation 
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length Lq, taking into account relativistic length contraction of a moving string in a box of 
size Lq! 

P = —F3— = 9j2 , (17) 

where g = (1 — f ^)~^/^ is the Lorentz factor associated with the motion of the string. We 
can then convert changes in string energy density to string accelerations via the equation 

dp dv V , ^ 

= d-vTt = jT^r ■ 

We have one such term for each string of type 1 and 2 and we require that they balance the 
energy gain due to the production of bridges, that is 

Pl.acc + P2,acc = Pl,2^3 ■ (19) 

Since the strings can have different tensions, the weighting depends on the relative string 



tensions |83|, so we introduce two weighting factors Wi{fii, fi2), ^2(^1, ^2) such that 
Wi{fii, 112) + u'2(/ii, 1^2) = 1 and require 

Pl,acc = Wi(/ii,/i2)pi,2^3 (20) 
P2,acc = W2{fil,fl2)Pl,2^3 ■ (21) 

Using Eqs. ([H]) and ([H]) we find 

Vi = d{l-Vi)wi{fii,fi2) — — ^jY^ 2 = l,2andi7^j (22) 

which are the new terms we seek. The velocity evolution equations for our model are 
therefore 

2-, f ki a ~ -.v fi3 £(t) \ 

Vi = {1 - Vi) { — - 2-Vi - dwi{fii, fX2) To-h I = 1,2 audi (23) 

\Ri a Vi fii L| J 

- - (I - ■ (2^) 

This picture is not complete, as we have not specified what happens when a string of 
type 3 collides with strings of type 1 or 2. For this toy model we will simply assume that 
the strings simply pass through one another in such an encounter, but more complicated 
models of this type will be considered in section IIVI Our aim here is to describe how such 
multi-string models may be constructed and to demonstrate that the presence of different 
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FIG. 4: Evolution of 7 = L/t and v for the three components of the network (bridge model). Type 
1 strings are shown in black, 2 in red and 3 in dotted blue. Type 3 are more heavily populated 
as they are continuously produced by interactions. Strings of type 2 are the less populated, being 
more massive than type 1, fj,2 = 2/ii. 



string types in the network, interacting with each other through terms hke flMj) . does not 
necessarily obstruct string scahng. Indeed, taking all loop production parameters Ci, C2 
and C3 equal to the canonical value c = 0.23 (coming from Nambu-Goto simulation of a 



single string network 



531 ]) and solving Eqs. f|T5|) - f|T6|) . fl23|) - fl24l) numerically, we find scaling 



behaviour for all three string types. This is illustrated in Fig.|l]for /i2 = 2/ii, /is = y/if + /il- 
For the weighting factors we have taken wi = /i2/(Aii + Ai2) and W2 = + 7x2), since 

the light string's motion is affected more due to its smaller inertia. We see that in this 
model, type 3 strings (shown in a blue dotted line) have smaller 7 and therefore become 
more abundant than types 1 and 2. This is simply a consequence of the fact that we have 
not allowed the breaking of type 3 strings into types 1 and 2: type 3 strings are continuously 
produced by interactions between type 1 and 2 strings, while there is no such mechanism 
to produce string types 1 or 2. Also note the important differentiation between string types 
1 and 2. This is because of the velocity dependence (see Eq. fl2^ ). since type 1 strings 
with less inertia lose more kinetic energy in an interaction. Lower velocities imply lower 
interaction rates for loop production which enhances the overall density of the light string. 
From the point of view of string density evolution, a scaling solution is generically reached. 
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as can be seen directly from the corresponding 73 = L^/t equation: 



73 ~ 2t 



73 7i72 t 



(25) 



Each term in the right hand side of fl25|) has a different dependence on 73. Further, the loop 
production term tends to increase 73, but as this happens this term weakens. Similarly, the 
term corresponding to the production of bridges tends to reduce 73 and in doing so it also 
weakens. Thus, there is always a value of 73 for which the right hand side of (!25|) is zero 
(scaling) and the dynamics is such that the system is driven towards that value. However, if 
we allow interactions between type 3 and type 1 or 2 strings by including bridge production 
terms in equations (fT5|) . then decreasing of 71, 72 in the bridge production term of equation 
f l25|l could in principle lead to a frustrated network. It is therefore important to consider 
more complicated situations separately, which we will do in section HVl 



3. Model with Zippers 

Now consider the case where collisions between strings of type 1 and 2 produce type 3 
segments through interactions of the 'zipper' type (Fig. [3]). In this case the original strings 
coalesce along their own length, so there will be extra terms in the corresponding energy 
density equations, describing the energy loss associated to the length of the colliding string 
which was converted to type 3. The energy density evolution equations thus become: 



Pi = -2-{l + Vi)pi FIT^' 1 = 1,2 (26) 



P3 = -2-{l + V^)p3 + r2T2 ■ (27) 



a L3 L1L2 

In general pi + P2 P3 so there will be an energy density of 



i(t) 

5pint = dv{pi + P2- /is) 7272"^^ (28) 



left from the interaction. We impose energy conservation by demanding that this energy 
excess is given as kinetic energy to the produced zippers, that is we require 

Ant = P3,acc • (29) 
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FIG. 5: Evolution of 7 = L/t and velocity v for the three components of the network (zipper 
model), for the case where the heavy (type 3) strings do not interact with the lighter ones. Type 
1 and 2 strings are shown with a black solid line, and type 3 is shown with a blue dotted line. All 
three components reach scaling with type 3 having smaller 7 (and so higher density) due to its 
continuous production from zipping of type 1 and 2 strings. 



This corresponds to the introduction of a new source term in the velocity equation of type 
3 strings. The velocity evolution equations for this model are then 



Ri 



a 
a , 



1,2 



, 2n (h ~v /xi + /i2 -/i3 Ll^jt) ' 



(30) 



(31) 



t3 u /is uii^2 

As explained above (section llll A ip . imposing energy conservation at zipping is a conserva- 
tive approach, which corresponds to switching off the possibility of damping energy through 
particle production or gravitational radiation. In this way we avoid the possibility of arti- 
ficially throwing away too much energy, which could lead to an unphysical scaling regime. 
The main purpose of this paper is to demonstrate that loop production alone can still pro- 
vide an efficient damping mechanism in the presence of links, so that we will assume energy 
conservation throughout the paper. When using phenomenological models of the type we 
present here to fit simulations, one can easily relax this assumption by introducing a tunable 
unbalance through reduction of the coefficient d in equation (13T|) . 

Solving Eqs. fl26l) - fl27j) . fl30|) - fl3T|) numerically with ci = 62 = £3 = 0.23 we find scaling 
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for all string types (Fig. \5^, where, as before, the heaviest type 3 strings are more highly 
populated. Again, this is because we have not included the possibility of 'unzipping' of type 
3 strings through their collisions with type 1 or 2: string length is continuously lost from 
type 1 and 2 strings into forming type 3 zippers, but the type 3 network only loses length 
through loop production. Allowing for the possibility that strings of type 3 collide with type 
1 (resp. type 2) producing type 2 (resp. type 1) segments one also finds scaling solutions 
were the most massive type 3 strings are less heavily populated (Fig. E]). We see therefore in 
this toy model that all string interactions have to be taken carefully into account, in order 
to obtain the correct scaling. We will construct more complicated models of this type in 
section |Vl 
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FIG. 6: An example of the evolution of 7 and v for a network (zipper model) where the heavy type 3 



strings can interact with the lighter ones and unzip. The heavy type 3 strings with jj^ = ^/xi + ^2 
are shown with a dotted blue line and the lighter strings with a black solid one. For simplicity we 
have taken /ii = ^2- All string types reach scaling but now type 3 strings can have smaller density 
due to their unzipping. 



B. General Case 

Arbitrarily complex models of multi-string networks can be built in a similar way by 
replicating the VOS equations ([S]),([n]) and adding extra terms of the type discussed above. 
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The result is a tower of ODE's of the form 



7-27-2 + 

a,k ^a-^i b,a<b 



^l^b 



(32) 



Vi 



bridges 



Via ^J^k i-ait) 



rLi a , Vi fli 



zippers 



5 Vab {^a + f^b- /ii) 



(33) 



b,a<b 

where = denotes the efficiency parameter for the process in which strings of type i and 
j interact to produce a type k segment, Vij is the average relative velocity between strings of 
type i and j, whereas ^fj(t) is the average length of links of type A;, produced by interactions 



between strings of types i and j around time t. The first sum in Eq. fl32p represents the 
energy lost from network i due to the length of type i strings that coalesced with other types 
to produce zippers. Note that no length from the colliding strings is lost when a bridge is 
produced (see section IIII A2|) so the sum is constrained only on interactions of the zipper 
type. On the other hand, the second sum models the energy gain of network i due to the 
production of both bridges and zippers of type i. Similarly, the first sum in Eq. (!33|) describes 
the slowing down of strings of type i due to their attachment to a newly formed bridge when 
they interact under the bridge configuration, while the second sum corresponds to the kinetic 
energy that must be given to zippers of type i to ensure energy conservation. More general 
string interactions can be modelled as combinations of the bridge and zipper types, and the 
relevant weighting can be encoded in the coefficients d^y Also, as mentioned above, one can 
relax the assumption of energy conservation at strin g in teractions, to include the possible 
effect of energy release through particle production [6^, [t^, by modifying some of these 
coefficients so that string mass-lengths and kinetic energies are not precisely balanced. 

Working as in sections III A H Hi A 21 we introduce the functions 7^ = Li/t and rewrite Eqs. 
( I32|) - (l33|) in the following form 



7i 



1 

2t 



CiVi 



zippers 



li 



ki 
li 



2/3vi 



a,k 

bridges 

- E 

a,k 



ia^ia o 
Ta 



1 



E d\bVab 



It 



ab 



it) 



b, a<b 



ilib 



(34) 



E7fc / \ Via f^k 1 ^ia (t) 
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lllb t 



(35) 
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Eqs. can be applied to any given model of multi-string network evolution, specific 

examples of which will be considered in sections IIVI and IVl 

Before moving to these applications, we comment on the size of segments i{t) produced 
in non-abelian interactions which is, in principle, another variable. One expects this to be 
a function of the correlation lengths of the colliding strings as well as the string velocities. 
In Ref. |40| . £{t) was expressed in terms of a constant 'zipping velocity' v^ip, as i{t) =fzipt. 
In Ref. j4l| on the other hand, all string types have the same correlation length L and 
all interactions are of the zipper type, so one can safely assume £{t) ~ L. Here, each 
string species has its own correlation length so we cannot adopt this approach. Instead 
we will consider bridge-type and zipper-type interactions separately and propose physically 
motivated ansatze. 

For bridge interactions, the original strings do not lose length but are slowed down due to 
the tension of the bridge formed between them. Therefore, given the tensions of all strings 
involved and the velocities and correlation lengths of the colliding segments, one can impose 
local energy and momentum conservation to find the length of the produced bridge. For 
comparable string tensions and correlation lengths, the following approximate formula holds 
for the length of the produced bridge: 

£b ^ ^{fiiL^vl + fi2L2vi) . (36) 

This is the formula we used in our bridge toy model, and we will also use it in section IIVI 
where we will study the evolution of Z^- strings with similar tensions, through a bridge- 
type model. Note that this is the maximum possible length that can be created because all 
the kinetic energy of the parent strings is lost in stretching the bridge to a length is- By 
integrating over different string orientations, it is clear that the average bridge length will 
be significantly less than in equation fl5^ . which corresponds to the parameter d'^j being 
significantly smaller than unity {d^j < 1). 

In the zipper case both of the interacting strings lose the same amount of length (which 
also equals the length of the produced zipper) and, since string direction changes after 
correlation length distances on the string, this cannot be larger than the smallest of the 
two correlation lengths. We could choose i{t) = min(Li,L2) but this would not be easy to 
implement in the equations. Instead, we will take 

n L1L2 



L1+L2 
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(37) 



a simple expression which returns a value smaller than, but not far off, the smallest of the 
two correlation lengths. This is a good approximation to the smallest correlation length if 
Li, L2 differ by an order of magnitude or more, and it returns half the correlation length 
for Li = L2. This is not however a problem, as the corresponding term in the evolution 
equations comes with a coefficient (i*^-, a free parameter of the model in which we can absorb 
this normalisation. The important issue is to make sure that C. has the correct scaling. We 
will also use this ansatz when modelling (p, g)-strings with zipping interactions in section IVl 
For more general interactions, interpolating between the two, the ansatz for ^{t) will 
depend on the specific model under study. Finally, for the collision velocity Vab between 
strings of type a and h we will take the magnitude of the relative velocity vector Vq — v^, 
averaged over directions, namely Vab = v'^a + "^t- 



IV. APPLICATION TO Zn STRINGS 

Consider a situation in which a continuous, simply connected group G is spontaneously 
broken to Z^q. Let 9 (where < < 2tx) parametrise a closed curve in physical space and 
denote the vacuum at position 9 by |^^). Then there is a group element g{9) G G which maps 
the vacuum at to that at 9, that is \9) = g{9)\0). Now consider the unbroken group Zn at 
9 = with elements hi = l,h2, ...,h]sf. The curve under consideration will encircle a string 
if g{2TT) = hi ^ 1. One would like to associate a different type of string to each non-trivial 
group element of the unbroken group. However, the group structure of Z^ imposes relations 
among various elements, which, under certain circumstances, allow one to identify strings 
corresponding to different group elements. 

To see this, consider the cyclic group Zn, which only has one generator h and elements 
hi = 1, h, h"^, h^~^. Acting on hN = h^~^ with the generator h gives back the identity, 
that is h-h^~^ = h^ = 1. By rewriting this as h^~^ = h~^ we see that strings corresponding 
to elements h2 = h and /iat = h^"^ have opposite fluxes. In some cases, depending on the 
details of the symmetry breaki ng p rocess, such strings and antistrings have been found to 
be connected by monopoles [tI, 72], and the evolution of such monopole-string networks has 
been discussed in Ref . jo^] . However, if the symmetry breaking process does not give rise to 



monopoles, each string and its corresponding antistring are topologically equivalent 



66 



72|. 



In particular, groups Z2 and Z^ give rise to one type of string only, while the lowest order 
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cyclic group which admits more than one distinct strings is Z^. 



Note however that branching can occur in Z 



networks too 



661], even when there is 



effectively only one type of string (see also Ref. 681). To see this, consider a volume element 
in physical space and assume that an /i2 string enters from one side. Charge conservation 
requires that the same amount of flux must leave the volume and this can be satisfied, 
for example, if an /i2 string also exits the volume from the other side, corresponding to 



the situation of a single string transversing the volume element [8J]. However, since 
(^2 = /ig = h!^ = h = h2, flux conservation is also satisfied if two /i2 strings enter from 
the other side, corresponding to three incoming /i2 strings joining at a vertex and forming a 
Y-type junction. But = so this can also be seen as an h2 string branching into two 
(outgoing) hs strings, two /i2 strings joining to an h^, or three hs strings emerging from a 
3-string vertex. From the point of view of string evolution [6'i1], /i2 and h-^ strings can be 
identified and the above configurations are equivalent. 

A. Modelling Z]\j networks 

In order to apply our string evolution model to a given string network, we need to study 
the topologically allowed outcomes from interactions between different types of string. This 
will allow us to set some of the d^j^. of equations (1M|) - (135|) . namely those corresponding to 
topologically forbidden interactions, to zero. Consider for example the case of Z^, which has 
elements l,h,h'^,h^. There are three types of string corresponding to /i2 = h, = h? and 

= h^. Since h ■ h = h? two strings of type 1 can join to a single type 2 string. Thus, 
two intersecting type 1 strings can either exchange partners or form a type 2 string segment 
(bridge) between them (Fig. [2]). Similarly, type 1 and type 2 strings can join to a single 
type 3 string, as h- h"^ = (Fig. [3]). Also, since h"^ -h? = 1, two type 2 strings can end on a 
vertex and so, depending on the details of the relevant symmetry breaking, either they are 



linked by a massive monopole or kink [t!], or they are self-conjugate |[72j]. In either case, the 
crossing of two type 2 strings can have two outcomes, as shown in Fig. [71 The presence of a 
massive monopole at the vertex where the strings meet does not directly affect macroscopic 
string evolution. Since monopole creation requires energy, such a configuration can result 
from high energy collisions only. This can be taken into account in our models by choosing 
an appropriately small efficiency parameter for such energetically disfavoured processes. As 
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FIG. 7: For self-conjugate strings, reconnection can happen in two distinct ways, as shown. 



with type 1 strings, two type 3 segments can join to a single type 2, because ■ = K^. 
Finally, /i ■ /i^ = 1 so a type 1 and a type 3 string can also end at a vertex. Then strings of 
types 1 and 3 can be seen as the same string, with opposite orientation. 

In the above example of strings we see that, although two colliding type 1 strings 
can form a bridge of type 2 string, intersecting type 2 strings can only exchange partners. 
Thus we have ^ 0, = in flM|) - fl35|) . Also, the intersection between strings of type 
1 and 2 can only lead to the formation of a type 1 bridge {d\2 = '^21 7^ 0) but not type 2 
{d\2 = (^21 — 0)- The parameters Cj and d^jp. can be grouped in a matrix 



A discussion of the topologically allowed outcomes of string collisions and the corresponding 
parameter matrices for the cases of Z^, Zq and Zj can be found in the appendix. 

For non-abelian string networks, string-string interactions are governed predominantly by 
topological, rather than energetic, considerations, except in extreme parameter regimes (e.g. 
strongly Type I strings). A link is produced between two colliding segments because this is 
the configuration which conserves charge, not because it minimises energy. Energy gains or 
losses associated to the colliding strings converting their length to the new string type are 
only marginal in general. Thus, one expects that such interactions can be approximated by 
those of the bridge type, in which the interacting strings are not losing length but the link 
forms and slows them down due its tension. We will now apply a bridge-type model to the 
above cases of Z^- strings. 




(38) 
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FIG. 8: Evolution of 7 and v for a Z3 network, where there is only one type of string, for d = 0, 
0.1c and 0.65c. The string network reaches scaling with smaller 7 as d increases. This corresponds 
to a larger string energy density as more Y-type junctions are being produced. 

1. Z3 strings 

The simplest case allowing the formation of Y-type junctions is Z^, where there is only one 
type of string and the collision between two string segments can either lead to reconnection 
or to the production of a bridge. The corresponding VOS model describing the evolution of 
such a network is given by the equations (see Eqs. fl34|) - fl35|) ) 



7 
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2t 



2p{l + v' 



2 + ^-4^ 
7 7^ t 
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(39) 



(40) 



where v = \/2v is the average relative velocity of strings and all other parameters/variables 
are as described in the text. Solving this system with c = 0.23 we find (Fig. [H]) scaling 
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solutions for d < 0.7c. The larger d is, the more entangled the network becomes and 
the scaling solution moves from its 'abelian' value to one with higher string density. For 
d ^ 0.7c the network becomes so dense that the bridge production term dominates, driving 
the system to an even denser state and thus spoiling scaling. In view of the discussion 
of last paragraph of IIII A2I this happens in the case of because both the interacting 
strings and the produced bridge are of the same type, so the dependence of the bridge 
production term on 7 changes from 7^ to 7"^. Then if this term gets to dominate, the string 
density increases (7 decreases) and the term becomes even stronger, so that scaling cannot 
be achieved. However, since intercommuting is energetically favourable compared to bridge 
production, one expects d < c. Then the formation of Y-type junctions increases the string 
density (Fig. [8j) but does not spoil the scaling property of the network. 
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FIG. 9: Evolution of 7 and v for the two components of a Z4 network, for ci = C2 = 0.23, 
dfi = O.OSci and = ci. Type 1 strings are represented by a solid black line and type 2 by a 
dotted blue one. We observe that type 1 strings are more abundant, since, unlike type 2, they can 
also be produced by interactions between type 1 and type 2 strings. 



2. Z4 strings 

For Z4 there are two types of string, say type 1 and type 2, and the corresponding 
parameter matrix is given by fl38p . A Z4 network can therefore be described by two copies 
of equations (!511) - (I5^ with non-zero parameters Ci, £2, dl^ and d\2 = dl^ (and with all 
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zipper terms set to zero). Again scaling solutions can be found when the bridge production 
terms do not dominate. For type 1 strings, since reconnection is energetically favourable 
compared to bridge production, we will assume dl^ < ci. For collisions between type 1 and 
2 strings, however, reconnection is not possible and the production of a type 1 bridge is the 
only option, so that can be large. One then finds scaling as long as does not exceed 
a critical value, which depends on the chosen d^i- For example choosing Ci = C2 = 0.23, 
df-^ = O.OSci (Fig. [H]) we find scaling, for both string types, for d\2 ~ 1.5ci. 



3. strings 



networks also consist of two types of string but with different interaction rules (see 
appendix). The parameter matrix is given by flA4p and the relevant VOS model arises by 
the appropriate adaptation of two copies of equations flM|l - fl35|) . Again, both string types 
achieve scaling when dfj are not larger than c (Fig. [TOj) . 
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FIG. 10: Evolution of 7 and v for the two components of a Z5 network. In this example we 
have chosen ci = 0.23, C2 = 0.3, df^ = 0.05ci, ^22 = 0.05c2, = 0.6ci and di2 = 0.4ci. 
Choosing ci = C2 and d\2 = £^12 instead leads to identical evolutions for the two string types, as 
the corresponding evolution equations are symmetric. 
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4- strings 



In the case of Zq there are three distinct types of strings (see appendix) with parameter 
matrix (]A6|) . The system is modelled by three copies of (lMI) - fl35l) . All three string types 
reach scaling for small enough fife's (Fig. [TTl) . 
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FIG. 11: Evolution of 7 and v for the three components of a Z6 network. Here, we have chosen 



ci = C2 = C3 = 0.23, = = O.OSci and = di2 = c^ia = O.lci. Type 1 strings are 



represented by a sohd black line, type 2 by a solid red, and type 3 by a dotted blue line. 



5. Zj strings 



Here there are also three string types with parameter matrix ( IA8|) . The relevant scaling 
solutions are shown in Fig. [121 



V. THE CASE OF COSMIC SUPERSTRINGS 



In this section we apply our model to cosmic superstrings, typically produced at the 



end of brane inflation 



3iv our moae 

a Q, Q, Q 



In this picture, a brane- ant ibrane pair are moving 



towards each other under their attractive interaction, and in doing so they give rise to an 
inflationary phase. Inflation ends when the branes collide and annihilate, leaving behind a 



network of D- and F-strings 
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29|. String interactions can lead to the formation 
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FIG. 12: Evolution of 7 and v for the three components of a Z7 network. We have chosen 
ci = C2 = C3 = 0.23, Jfi = 0.05ci and = ^12 = ^13 = ^13 = ^^22 = ^3 = ^23 = ^3 = 0-lci. As 
before, type 1 strings are represented by a black sohd Hne, type 2 by a red solid and type 3 by a 
blue dotted line. 



of bound states between p F-strings and q D-strings, referred to as {p, g)-strings, where p, q 
are coprime |851]. This situation corresponds more closely to interactions of the zipper type, 
rather than bridge production, as the strings coalesce along their own length in forming 
a bound state. Also, unlike strings, where different string types could have the same 
tension, here each type of string has different tension. F-strings, being perturbative objects 
have tensions proportional to the square root of the string coupling, namely 



(41) 



while D-strings are non-perturbative and have 



(42) 



In flat space the tension of {p, q) strings is given by 



(43) 
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A. Modelling superstring networks 



To model cosmic superstring networks we consider two types of elementary string, type 
1 with tension /ii = jiQ^fgl (corresponding to the F-string) and type 2, tension /i2 = /io/y^ 
(the D-string). When two elementary strings of the same type collide, they may pass through 
one another or reconnect by exchange of partners. However, a type 1 and a type 2 string 
can bind together to form a (1, ±1) string segment. In general, when a (p, q) and ajp'^o') 



string collide they can bind in two ways, depending on the angle of collision [2c 
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4l|, 



forming either a {p + p' ,q + q') or a (p — p', q — q') segment (Fig. [T3|) . where without loss of 
generaUty we have taken p > p'. G.ven that btnding has taken place, the ptobabUity that 

the additive/subtractive process has occurred is given by [41j] 

pp'gl + qq' W 



(44) 



(p2^2 + ^2)1/2 (p.2^2 + ^,2)l/2yy' 

where we have assumed that the RR scalar is zero. Assuming that the zipper type interaction 
is a good approximation to this binding, we can model such a network by equations flMl) - fl5^ 
with zero bridge terms and 



{p±p',q±q') 

{p,q),{p',q') 



with '^(p,q),(p',q')'s range found in Ref. [29|. 
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FIG. 13: A {p,q) and a {p',q') strings, can bind together in two distinct ways, forming either a 
(P + P\q + q') or a (p - p', q - q') zipper. 
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Fig. [13] shows the scahng density results of a model containing bound string states up 
to (3, ±1) and (1,±3). In agreement with Ref. |41| the string densities of heavier states 
fall rapidly as the string tension increases. This is because of equation (j331), which tends 
to give a small value of (hence a large value of P^) for large p's and g's, so that 
heavy strings have the tendency to break to lighter ones. It is therefore possible to obtain 
a very good approximation of such a network by truncating the equations at a relatively 
low s = p + g, like, in this case, s = 4. As one reduces the intercommuting probability 
(Fig. [H]) string densities increase, as expected, and the fall of density with increasing string 
tension becomes more prominent. The relative importance of light to heavier string states 
depends on both the string coupling Qs (appearing in the string tensions and -P^), and the 
intercommuting probability absorbed in c's and (i's. However, the generic behaviour is that 
the string number density is dominated by strings of type (1, 0), (0, 1) and (1, ±1), all having 
comparable number densities, while higher composite states are heavily suppressed. 
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FIG. 14: Evolution of normalised string density N = pt^ / ji for the lightest components of a cosmic 
superstring network. F-strings are shown with a black solid line, D-strings with dotted blue, (1, ±1) 
strings are in green, (2, ±1) in red and (3, ±1) in yellow. Type (1, ±2) and (1, ±3) are shown in 
dotted red and dotted yellow respectively. All string components are seen to approach scaling with 
heavier bound states less abundant. As one reduces the intercommuting probability (left to right) 
the string densities increase and light strings become more dominant. 
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The fact that the number density of the first bound state (1, ±1) is comparable to that of 
the unbound strings can be understood in terms of the difference in string tension between 
F- and D-strings. This makes it energetically favourable for the light (1,0) strings to be in a 
bound (1, ±1) state. One can then envisage such a network as being predominantly composed 
by a heavy component of D-strings, together with a 'cobweb' of light F-strings, which tend 
to stick to the D-string network, giving rise to another heavy FD-string component. 

In our velocity dependent model we separately evolve the rms velocities of each string 
component of the network, which, since each string type has a different correlation length, 
can be significantly different. This is shown in Fig. [15] where, evidently, the less abundant 
strings also have smaller velocities. This could in principle be observable as the velocity 



enters the temperature discontinuity through the Kaiser-Stebbins effect 
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FIG. 15: Evolution of the rms string velocity v for the first lightest components of a cosmic 
superstring network. F-strings are shown with a black solid line, D-strings with dotted blue, 
(1, ±1) strings are in green, (2, ±1) in red and (3, ±1) in yellow. The dotted red and yellow curves 
correspond to strings of type (l,ib2) and (l,ib3). All string components asymptotically reach a 
constant value, which increases as one reduces the intercommuting probability (left to right). 

The above results provide further evidence in favour of the possibility of scaling in cosmic 
superstring networks and other 'entangled' networks with junctions. Particularly interesting 
is the fact that the scaling solutions we found do not rely on any additional energy damping 
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mechanism, apart from loop production through self-intersections of strings. Indeed, by im- 
posing energy conservation during zipping events, we have ensured that no energy is damped 
during the zipping processes, which would facilitate the reaching of a scaling regime. Pos- 
sible non-elasticity of the zipping events would provide an additional decay channel, which 
would therefore further support our scaling results. One should keep in mind that the mod- 
els we presented here and applied to the case of cosmic superstrings, are phenomenological 
in nature and, as such, do not model in detail a number of microphysical processes, related, 
for example, to the dynamics of string junctions ^4] (see Ref. 75] for recent field theory 
simulations). However, progress in understanding the complex, non-linear evolution of cos- 
mic strings has been possible with the combination of different approaches (microphysical 
models, Nambu-Goto and field theory simulations, phenomenological models), and there is 
already strong evidence that models like those presented here capture most of the relevant 
physics. Being able to construct such analytic models of multi-tension string evolution and 
come up with quantitative predictions is an important step in our understanding of non- 
abelian cosmic (super) string evolution. It would be interesting to combine these models 
with other approaches to cosmic superstring modelling as for example field theory models 
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461] or Nambu-Goto multi-tension simulations. 



VI. CONCLUSION 



We have presented analytic velocity-dependent models for the cosmological evolution of 
non-abelian string networks. Apart from ordinary abelian intercommutations, these models 
account for interactions between different string types, producing Y-type junctions with 
linking segments stretching between the originally colliding strings. We have described such 
interactions in terms of two limiting cases, namely zipper-type interactions, where the link 
is produced by the zipping of the colliding strings into a bound state, and bridge-type ones, 
where the link is a newly formed string of a third type and the parent strings do not lose 
significant string length. A general Y-junction forming interaction can be describing by 
a combination of these two inter-related mechanisms, and the relative weighting depends 
on the type of strings under study. In non-abelian field theory, where the interactions are 
topologically constrained, one can argue that the bridge-type picture is a good approximation 
as the colliding strings do not gain energy by converting length to the new type. On the 
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other hand the zipper interaction corresponds closely to the case of cosmic superstrings, 
where D-strings and F-strings tend to bind together to form heavier composites. 

In the bridge case, we have applied our models to the evolution of Zn string networks and 
found that the presence of Y-type junctions does not generally lead to string frustration. 
Instead, scaling solutions exist for a wide range of physically relevant choices of parameters. 
Here, we have thus discussed the physical properties and scaling behaviour of Z]\r networks 
for = 3, 4, 7 noting the same qualitative behaviour. 

In the zipper case, we have modelled a cosmic superstring network as one consisting of 
two types of elementary strings (corresponding to F- and D-strings) which can zip together 
to form bound states between p strings of one type and q strings of the other. We have 
demonstrated scaling of all st ring types, with heavier strings generally less populated than 



lighter ones, as noted in Ref. 



4l| . We have also obtained the scaling velocities of each net- 



work component, with heavier strings moving slower than lighter ones. The general picture 
emerging from this analysis, is a string network whose number density is dominated by the 
lightest (unbound) strings and the first bound state between them, with all heavier bound 
states being suppressed. The first bound state develops a comparable number density to 
the unbound strings, even though it has a higher string tension. This is to be understood 
in terms of the difference in string tension between F- and D-strings, which makes it ener- 
getically favourable for the lightest F-strings to be in a bound FD-state. Thus, a scaling 
superstring network is dominated by D-strings and the first FD bound state, with a cobweb 
of lighter F-strings of comparable number density, but subdominant energy density. 

The strength of our scaling results stems from the fact that we have not allowed additional 
energy damping processes associated to the production of junctions (though our models 
can accomodate these possibilities also), as we have enforced energy conservation in the 
relevant string interaction events. Thus, first, we can be sure that we are not observing a 
'spurious' scaling behaviour due to artificially throwing away some of the energy involved, 
and, second, our scaling results provide evidence that loop production alone may be sufficient 
for scaling (even in the presence of junctions) without the necessity of extra energy losses, as 
for example massive radiation at zipping. However, it should be highlighted that our findings 
provide evidence, rather than proof, of scaling in such networks. Indeed, our models are 
phenomenological in nature and, as such, have their own limitations. In particular, the 
detailed microphysics of junctions, for example the possibility of unzipping, is not modelled 
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explicitly, but only statistically and indirectly, through choices of the relevant parameters. In 
principle, if juncion dynamics is sufficiently complicated, these coefficients could inherit time 
dependence, a possibility we have not considered. One should therefore combine our results 
with other complementary approaches to this difficult problem, each of which captures only 
part of the relevant physics. What is encouraging, is that all complementary approaches 
seem to independently point towards scaling in these networks. 
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46| 



It would be interesting to compare these models to field theory simulations 
and see whether quantitative agreement can be established. It would also be desirable to 
generalise them to include a second length-scale for each network component, as, for small 
intercommuting probabilities likely to appear in cosmic superstring networks, a two-scale 
VOS model (see 



Goto simulations of single string networks 



76| . Other variants of our models could also 



be constructed. For example, the models presented here assume that the production of 
junctions can be approximated by the two mechanisms described in Ref. 3l(|, namely the 
'zipper', where the string lengths lost from each of the zipping strings equal the length of 
the produced link, and the 'bridge', where the lengths of the colliding strings are preserved. 



Recent studies of junction formation based on the Nambu-Goto action [43|, |4J] suggest a 
slightly different mechanism where strings lose or gain length subject to a constraint, which 
the vertex has to satisfy. This mechanism could be readily accommodated in our models. It 
is an interesting question to what extent scaling results depend on such choices. Analytic 
models like these will be useful in comparing and contrasting with other approaches to 
'entangled' network evolution. Although quantitative agreement with numerical simulations 
has not yet been investigated, these models provide further evidence that scaling behaviour 
in networks with junctions is possible, if not generic. 
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APPENDIX A: PARAMETER MATRICES FOR Zn STRINGS, 4 < AT < 7 



In this appendix we explore the topologically allowed outcomes from string collisions in 
the case of Z4, Z5, Zq and networks. 

1. Z4 strings 

Consider Z4 with generator h and elements hi = l,h,h'^,h^. There are three types of 
string (say type 1, 2 and 3) corresponding to the non-trivial group elements h, h"^ and 
respectively. We have h ■ h — h? so that strings of type 1 can produce a type 2 segment 
when they interact, that is 

1 + 1^2. (Al) 
Thus, the crossing of two type 1 strings can have two different outcomes, as shown: 




Similarly, we have: /i-/i^ = /i^^l-|-2— >3, corresponding to 




• = 1 =r- 2 + 2 — > 0, meaning that two incoming type 2 strings can end on a vertex 
and so the interaction between type 2 strings can then have two outcomes: 




/i^ • /i^ = /i^ 3 + 3 — > 2, which is as in the first diagram, 
and: 
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= 1^1 + 3^0, that is a type 1 and and a type 3 strings can end on a vertex, 
and can be treated as being the same string type with opposite orientation: 





Hence there are only two distinct types of string, type 1 = —3 and 2, with type 2 being 
self-conjugate. The parameter matrix corresponding to the above interactions is (see also 
section IIV Al) : 



M 



^{~c,Jl,) (4,0) 

^ (4,0) (0,C2) 



(A2) 



2. Z5 strings 

For Z5 we have hi = l,h,h'^,h^,h^. Associating a string type to each non-trivial generator, 
as before, we have: 

1 + 1^2 
1 + 2^3 
1 + 3-^4 

1 + 4-^0^1^-4 1 + 1 

2 + 2-^4 1 + 2 

2 + 3^0 ^2 = -3 l + (-2) 
2+4^1 2+2 

3 + 3^1 

3 + 4^2 

4 + 4-^3 

Again, there are two types of string with interaction diagrams: 
1 + 1 -> 2 

^ 1 



-2 
-1 
-1 




(A3) 
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These are described by the parameter matrix: 



/ 



M 



{d2i,d2i) ('^22)^2) 



(A4) 



Note that there is no self-conjugate string in this case. 
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3. strings 

For we have hi = 1, h, h"^, h^, h'^, h^. Now: 

1 + 1^2 
1 + 2^3 
1 + 3^4 
1 + 4^5 

1 + 5^0 =^ 1 = -5 

2 + 2^4 
2 + 3^5 

2 + 4^0 ^2 = -4 

2 + 5^1 

3 + 3^0 
3 + 4 ^ 1 

3 + 5^2 

4 + 4^2 

4 + 5^3 

5 + 5^4 

Some of the new diagrams arising are 
1 + 2^3 
1 + (-2)^-1 



3 self —conjugate 



1 + 1 
1 + 2 

1 + 3 

1 + (-2) 

2 + 2 
2 + 3 

2 + (-l) 
3 + (-2) 

3 + (-l) 




2 
3 

-2 
-1 
-2 
-1 
-1 

1 

2 



(A5) 



and: 
2 + 2 
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There are only three string types with parameter matrix: 



I (gi,4,0) (4,0,4) (o,4>o)\ 

(4,0,4) (0,|,0) (4,0,0) 

"22 

\, (0,4,0) (4,0,0) (0,0, C3) ) 



String 3 is self-conjugate. 



4. Z-j strings 



For Z-j we have hi — 1, /i, /i^, /i^, /i^, /i^, /i^. Working as before: 



1 + 1 
1 + 2 
1 + 3 
1 + 4 
1 + 5 

1 + 6 

2 + 2 
2 + 3 
2 + 4 
2 + 5 

2 + 6 

3 + 3 
3 + 4 
3 + 5 

3 + 6 

4 + 4 
4 + 5 

4 + 6 

5 + 5 
5 + 6 



2 
3 
4 
5 
6 

4 
5 
6 

1 
6 

1 
2 
1 
2 
3 
3 
4 



^ 1 = -6 



3 = -4 



1 + 1 
1 + 2 

1 + (-2) 

1 + 3 

1 + (-3) 

2 + 2 

2 + 3 

2 + (-3) 

3 + 3 

3 + (-2) 

3 + (-l) 



2 
3 



1 

2 
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Again, there are only three string types, but no self-conjugate strings. The parameter matrix 
is now: 



M 



' (ci,dl,0) (^1,0,^1) (0,4,4)^ 
(Jl2,0,dl2) (0,C2,4) (4,4,0) 

v(o,4,4) (4,4,0) (4,o,c3) ] 



(A8) 
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